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Abstract. 

A model for opinion dynamics (Model I) has been recently introduced in which the binary 
O ' opinions of the individuals are determined according to the size of their neighboring domains 

O ' (population having the same opinion). The coarsening dynamics of the equivalent Ising model 

shows power law behavior and has been found to belong to a new universality class with the 
^ • dynamic exponent z = 1.0 ± 0.01 and persistence exponent 6 ~ 0.235 in one dimension. The 

, ^ ' critical behavior has been found to be robust for a large variety of annealed disorder that has 

^ . been studied. Further, by mapping Model I to a system of random walkers in one dimension 

with a tendency to walk towards their nearest neighbour with probability e, we find that for 
any e > 0.5, the Model I dynamical behaviour is prevalent at long times. 



1. Introduction 



> 

\^ ' Sociophysics has emerged as one of the important areas of research during recent times. The 

Tlj" . concepts of statistical physics find application to many situations that occur in a social system 

| with the assumptiou that individual free will or feelings do not take crucial role in these situations 

CD ■ [HE]- One of the major issues that has attracted lots of attention is how opinions evolve in a 

social system. Starting from random initial opinions, dynamics often leads to a consensus which 
means a major fraction of the population support a certain cause, for example a motion or a 
^ candidate in an election etc. 

^ I Simulating human behaviour by models effectively implies quantifying the outcome of the 

^ ■ behavior by suitable variables having continuous or discrete values. Different dynamical rules 

- - - are proposed for the evolution of these variables, depending on how these variables change with 

time following social interactions. Thus, a social system can be treated like a physical system. 
For example, in case of opinion dynamics, if the opinions have only discrete binary values, the 
social system can be regarded as a magnetic system of Ising spins. 

In this context, Schelling model [3j, proposed in 1971, seems to be the very first model of 
opinion dynamics. Since then, a number of models describing the formation of opinions in a 
social system have been proposed [4]. While on one hand these models attempt an understanding 
of how a society behaves and social viewpoints evolve, on the other hand, these provide rich 
complex dynamical physical systems suitable for theoretical studies. 

Dynamics of complex systems has become a subject of extensive research from several aspects. 
For many such systems, e.g., traffic or agent based models, one cannot define a conventional 



Hamiltonian or energy function. The only method by which one can study the steady state 
behaviour of such systems is by looking at the long time dynamics. Nonequilibrium dynamics 
involves the evolution of a system from a completely random initial configuration and associated 
with this evolution are several phenomena of interest like domain growth or persistence that 
have been studied, for example, in spin systems. Since in many sociophysics model, one can 
have variables analogous to spin variables, these phenomena can be readily studied here. An 
important objective is to identify dynamical universality classes by estimating the relevant 
dynamical exponents. 

Another point of interest in studying dynamics is that many systems may have identical 
equilibrium behaviour but behave differently as far as dynamics is concerned. For example, 
Ising spin dynamics with or without conservation belong to different dynamic universality class 
although their equilibrium behaviour is identical. 

Apart from the dynamical behaviour, different kinds of phase transitions have also been 
observed in these models by introducing suitable parameters. One such phase transition can be 
from a homogeneous society where everyone has the same opinion to a heterogeneous one with 
mixed opinions [5j. 

Change in the opinion of an individual takes place in different ways in different models. For 
example in the Voter model [6], an individual simply follows the opinion of a randomly chosen 
neighbour while in the Sznajd model [7J, the opinion of one or more individuals are changed 
following more complicated rules. 

In this article, we review the dynamical studies in a recently introduced model [8J (to be 
referred to as Model I henceforth) and its variants in which a new rule of updating (discussed 
in detail in section 2) is introduced. In one dimension. Model I can be visualized as an Ising 
spin chain (the binary opinions are represented by Ising spins) with a Glauber-like dynamics, 
where a spin only at the domain boundary can flip. In this new model, the state of the spin is 
determined by the state of the neighboring domain larger in size (detailed description is given 
in the next section). Model I shows strikingly different dynamical behavior compared to known 
models of opinion dynamics or spin dynamics. Several observables show power law decay and 
the exponents strongly suggest a new universality class. 

The introduction of disorder in various forms have also been considered which shows that the 
Model I dynamical behaviour is not affected by annealed disorder. We also report (in section 
3) new results when a mapping of this model is made to an equivalent reaction diffusion model 
where the walkers move towards their nearest neighbour with probability e. Once again Model I 
dynamical behaviour is seen to exist at long times for any e > 0.5 (which means a bias towards 
the nearest neighbour) showing the extreme robustness of the model. 

2. Description of the Model I 

In a model of opinion dynamics, the key feature is the interaction of the individuals. Usually, in 
all the models, it is assumed that an individual is influenced by its nearest neighbours. Model I 
is a one dimensional model of binary opinion in which the dynamics is dependent on the size of 
the neighbouring domains as well. Here an individual changes his/her opinion in two situations: 
first when the two neighbouring domains have opposite polarity, and in this case the individual 
simply follows the opinion of the neighbouring domain with the larger size. This case may arise 
only when the individual is at the boundary of the two domains. An individual also changes 
his/her opinion when both the neighbouring domains have an opinion which opposes his/her 
original opinion, i.e., the individual is sandwiched between two domains of same polarity. It may 
be noted that for the second case, the size of the neighbouring domains is irrelevant. When the 
two neighbouring domains are of the same size but have opposite polarity, the individual will 
change his/her orientation with fifty percent probability. 

The binary opinions can be represented by a system of Ising spins where the up and down 



states correspond to the two possible opinions. The two rules followed in the dynamical evolution 
in the equivalent spin model are shown schematically in Fig. [1] as case I and II. In the first case 
the spins at the boundary between two domains will choose the state of the left side domain (as 
it is larger in size). For the second case the down spin flanked by two neighbouring up spins will 
flip. 
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Figure 1. Dynamical rules for Model I: in both cases the encircled spins may change state; 
in case I, the boundary spins will follow the opinion of the left domain of up spins which will 
grow. For case II, the down spin between the two up spins will flip irrespective of the size of the 
neigbouring domains. 

The main idea in Model I is that the size of a domain represents a quantity analogous to 
'social pressure' which is expected to be proportional to the number of people supporting a 
cause. An individual, sitting at the domain boundary, is most exposed to the competition 
between opposing pressures and gives in to the larger one. This is what happens in case I shown 
in FiglU The interaction in case II on the other hand implies that it is difficult to stick to one's 
opinion if the entire (immediate) neighbourhood opposes it. 

Defining the dynamics in this way, one immediately notices that case II corresponds to what 
would happen for spins in a nearest neighbour Ferromagnetic Ising Model (FIM) in which the 
dynamics at zero temperature is simply an energy minimization scheme. However, the boundary 
spin in the FIM behaves differently in case I; it may or may not flip as the energy remains same. 
In the present model, the dynamics is deterministic even for the boundary spins (barring the 
rare instance when the two neighbourhoods have the same size in which case the spin flips with 
fifty percent probability). 

In this model, the important condition of changing one's opinion is the size of the neighbouring 
domains which is not fixed either in time or space. This is the unique feature of this model. In 
the most familiar models of opinion dynamics like the Sznajd model [7] or the voter model [6j, 
one takes the effect of nearest neighbours within a given radius and even in the case of models 
defined on networks [QJ , the influencing neighbours may be nonlocal but always fixed in identity. 

In the equivalent spin model, if -L_|_ is the number of up spins and L_ is the number down 
spins, the order parameter is defined as m = |L+ — L_|/L. This is identical to the (absolute 
value of) magnetization. Starting from a random initial configuration, the dynamics in Model I 
showed that it leads to a final state with m = 1, i.e. a homogeneous state where all spins have 
the same value (either +1 or -1). It is not difficult to understand this result; in absence of any 
fluctuation, the dominating neighbourhood (domain) simply grows in size ultimately spanning 
the entire system. 

In the spin picture, the dynamics can be described in terms of the movement of the domain 
walls and as the dynamics progresses, number of domain walls goes on decreasing. Monte Carlo 
simulations showed that the domain dynamics and the dynamics of the order parameter obey 
conventional power law variations: the fraction of domain walls D{t) oc t~^^^ with z = l.OOibO.Ol 
and order parameter m{t) oc i^/^^ (Fig. [5]). 
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Figure 2. Variation of the order parameter m with time for two different system sizes along 
with a straight Une (slope 0.51) shown in a log-log plot. Inset shows the decay of fraction of 
domain wall D with time. 



The persistence measure (i.e., probability that a spin has not flipped till time t) [TO] showed 
the familiar power law behaviour: P{t) oc t~^, where 9 is the persistence exponent. For finite 
system of size L, P{t, L) is known to behave as |1H [T2] 



and at large times, the persistence probability saturates at a value oc L~". Therefore, for x << 1, 
/(x) oc x~" with a = z9. For large x, /(x) is a constant. Estimate for both z and 9 using the 
above scaling form showed that 9 = 0.235 ± 0.003, and a z value 1.04 it 0.01 (Fig. [3|). 

Model I therefore showed a novel dynamical behaviour with values of z and 9 quite different 
from those of the one dimensional Ising model [13j and other opinion/voter dynamics models 
[m [151 [16] . Specifically in the Ising model, z = 2 and 9 = 0.375 and for the Sznajd model the 
persistence exponent is equal to that of the Ising model. 



P(t,L) oc t~^f{L/t 
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Figure 3. The collapse of scaled persistence probability versus scaled time using 9 = 0.235 and 
z = 1.04 is shown for different system sizes. Inset shows the unsealed data. 



2.1. Model I with a cutoff 

As a variation of the Model I, one can introduce a cutoff [17] denoted by a parameter p, such that 
the maximum size of the neighbouring domains an individual can sense is given hy R = pL/2 
in a one dimensional system of L agents with periodic boundary condition, p = 1 corresponds 
to the original Model I. A finite cutoff (i.e., p < 1) puts a restriction on the domain sizes which 
may correspond to geographical, political, cultural boundaries etc. The case with uniform cutoff 
signifies that all the individuals have same kind of restriction. The case with random cutoffs, 
which is perhaps closer to reality, has also been studied briefly. Once again one can represent 
the system by Ising spins. 

Whenever R is kept finite, that is, R does not scale with system sizes (implying p — )• in the 
thermodynamic limit), the dynamics leads to the equilibrium configuration of all spins up/down 
and the dynamic exponents also turn out to be identical to those corresponding to the nearest 
neighbour Ising values (i.e., = 0.375 and z = 2). 

The results with a nonzero p showed that there is a crossover phenomena in dynamics. For 
an initial time ti = pL/2, the dynamics remains Model I-like as the cutoff does not affect the 
dynamics till this time. Beyond this time, the dynamics becomes diffusive. Such a crossover 
is, however, not apparent from the log-log plots of the relevant dynamical quantities (Fig. H]). 
This is because the crossover occurs between two different types of phenomena. The first one 
is pure coarsening in which domain walls prefer to move towards their nearest neighbours as in 
Model I and one gets the expected power-law behaviour. There is a correlation length which is 
growing; as long as this length is smaller than R one will observe ballistic growth (as in Model 
I). When the correlation length is larger than R one observes domain diffusion. The diffusive 
behavior is not evident in the coarsening process because ballistic dynamics "leaves" the system 
into a non- typical configuration which is evidently far from those on diffusion paths. At ti, 
some special configurations are generated in which only a few domain walls remain, the density 
of domain walls going to zero in the thermodynamic limit. Therefore beyond ti, the second 
phenomenon involves pure diffusion of these domain walls which remain non-interacting up to 
large times. Only when the distance between two domain walls become < R the Model I like 
dynamics prevails again before the two are annihilated. This takes place over a much shorter 
time scale such that the physical quantities attain saturation values (m = 1 and D = 0) very 
fast in the last part of the dynamics (Fig. H]). 




Figure 4. Variation of the (a) order parameter and the (b) fraction of domain walls with time 
for different finite cutoff p = 0.1,0.2,0.4 (L = 3000). 



The only dynamic exponent in the diffusive regime is the diffusion exponent z = 2, which 
is distinct from the growth exponent z = 1. So the two dynamic exponents not only differ in 
magnitude, but they are also connected to distinct phenomena. Using both numerical data and 
some analytical arguments, it could be shown that the behaviour at t > ti is indeed diffusive. 
The saturation time to reach equilibrium was found to be 

tsat = apL + b{l-pfL^, (2) 

which also showed that for L — ?■ oo, z = 2. 

The persistence behaviour, however, does not show any power law behaviour corresponding 
to the diffusive behaviour, i.e., it does not decay like t~C).375 ^ usual reaction diffusion 

system. This is because of the special configurations in which the system is left at time ti. 
However, the long time behaviour showed that the persistence probability saturates in time and 
decays with the system size as with a = 0.235, which is the Model I value. In fact one 
can obtain a collapse for the persistence data using a = 0.235 in both the regions t < ti (with 
z = 1) and t > ti (with z = 2) for any p ^ 0. Thus the crossover phenomenon occurs with a 
novel characteristic behaviour. 

2. 2. Effect of disorder 

The Model I described so far has no disorder, which can be introduced in several possible ways. 
We discuss three cases in the following. 

2.2.1. Effect of rigidity parameter: Since every individual is not expected to succumb to social 
pressure, Model I can be modified by introducing a parameter p called rigidity coefficient which 
denotes the probability that people are completely rigid and never change their opinions [8j. 
This means there are pN rigid individuals (chosen randomly at time t = 0), who retain their 
initial state throughout the time evolution. Thus the disorder is quenched in nature. Such rigid 
individuals had been considered earlier in [18j. 

The limit p = 1 corresponds to the unrealistic noninteracting case when no time evolution 
takes place; p = 1 is in fact a trivial fixed point. For other values of p, the system evolves to a 
equilibrium state. 

The time evolution changes drastically in nature (compared to Model I) with the introduction 
of p. All the dynamical variables like order parameter, fraction of domain wall and persistence 
attain a saturation value at a rate which increases with p. Power law variation with time can 
only be observed for p < 0.01 with the exponent values same as those for p = 0. The saturation 
or equilibrium values on the other hand show the following behaviour: 

rris oc Af"°V~'^i 
oc p~^' 

Ps=a + hp-P^ (3) 

where in the last equation a is a constant ~ 0.06 independent of p. The values of the exponents 
are ai = 0.500 ± 0.002, pi = 0.513 ± 0.010, /32 = 0.96 ± 0.01 and = 0.430 ± 0.01. 

It can be naively assumed that the N p rigid individuals will dominantly appear at the domain 
boundaries such that in the first order approximation (for a fixed population), D oc 1/p. This 
would give m oc 1/y^ indicating (3i = 0.5 and /?2 = 1. The numerically obtained values are in 
fact quite close to these estimates. 

The results obtained above can be explained in the following way: with p 7^ 0, the domains 
cannot grow freely and domains with both kinds of opinions survive making the equilibrium 
rus less than unity. Thus the society becomes heterogeneous for any p > when people do 



not follow the same opinion any longer. The variation of with shows that nis — )• in 
the thermodynamic limit for p > 0. Thus not only does the society become heterogeneous at 
the onset of p, it goes to a completely disordered state analogous to the paramagnetic state in 
magnetic systems. Thus one may conclude that a phase transition from a ordered state with 
m = 1 to a disordered state (m = 0) takes place for p = . 

Since the role of p is similar to domain wall pinning, one can introduce a depinning probability 
factor p which in this system represents the probability for rigid individuals to become non-rigid 
during each Monte Carlo step, p relaxes the rigidity criterion in an annealed manner in the 
sense that the identity of the individuals who become non-rigid is not fixed (in time), p = 1, 
one gets back Model I whatever be the value of p, and therefore p = 1 signifies a line of (Model 
I) fixed points, where the dynamics leads the system to a homogeneous state. 
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Figure 5. The flow lines in the p — p plane: Any non-zero value of p with p = drives the 
system to the disordered fixed point p = 1- Any nonzero value of p drives it to the ordered state 
{p = 1, which is a line of fixed points) for all values of p. 

With the introduction of p, one has effectively a lesser fraction p' of rigid people in the society, 
where 

P=pil-p). (4) 

The difference from the previous model with quenched rigidity is, of course, that this effective 
fraction of rigid individuals is not fixed in identity (over time). Thus when p 0, p ^ 0, we 
have a system in which there are both quenched and annealed disorder. It is observed that for 
any nonzero value of p, the system once again evolves to a homogeneous state (m = 1) for all 
values of p. Moreover, the dynamic behaviour is also same as Model I with the exponent z and 
9 having identical values. This shows that the nature of randomness is crucial as one cannot 
simply replace a system with parameters {p ^ 0, p ^ 0} by one with only quenched randomness 
{p' ^ 0, p' = 0} as in the latter case one would end up with a heterogeneous society. We 
therefore conclude that the annealed disorder wins over the quenched disorder; p effectively 
drives the system to the /U = 1 fixed point for any value of p. This is shown schematically in a 
flow diagram (Fig. [5]). It is worth remarking that it looks very similar to the flow diagram of 
the one dimensional Ising model with nearest neighbour interactions in a longitudinal fleld and 
finite temperature. 

2.2.2. Effect of thermal-noise like disorder : In thermodynamical systems, the effect of thermal 
disorder is a highly important issue. In Model I which describes a social system, a disorder 
analogous to thermal noise has been introduced [19]. 



Let dup and ddown be the sizes of the two neighbouring domains of type up and down of 
a spin at the domain boundary (excluding itself). In Model I, probability P[up) that the 
said spin is up is 1 if dup > ddown, 0.5 if dup = ddown and zero otherwise. In the simplest 
possible way to introduce disorder, one may take the probability of a boundary spin to be up 
as P{up) = dup/ {dup + ddown)- However, there is no parameter controlling the stochasticity 
here and moreover, the results are identical to the original Model I and therefore this kind of 
stochasticity is not of much importance. 

In order to introduce a noise like parameter which can be tuned, the probability that a spin 
at the domain boundary is up was taken to be 

P{up) oc e^('^"f-'^''°™"), (5) 

and it is down with probability 

P{down) oc e^('''^°™"-'^''p). (6) 

The normalized probabilities are therefore P{up) = exp /3A/(exp(/3A) + exp(— /3A)) and 
P{down) = 1 - P{up), where A = {dup - ddown)- 
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Figure 6. Snapshots of the system in time for different values of /3 = 0.0, 0.005, 0.1, and /3 — >• oo 
(top to bottom) showing that even for very small non-zero /3, the system equilibrates towards a 
homogeneous state much faster compared to the /3 = case. These snapshots are for a L = 100 
system. 



Obviously, /3 ^ oo corresponds to Model I while letting /3 = we have equal probabilities of 
the up and down states, making it equivalent to the zero temperature dynamics of the nearest 
neighbour Ising model. Since the equilibrium states for the extreme values /3 — t- oo and /3 = 
are homogeneous (all up or all down states), it is expected that for all values of /3 they will 
remain so as is indeed the case. 



It is useful to show the snapshots of the evolution of the system over time for different (3 
(Figj6]): to be noted is the fact that for any non-zero /3 however small, the system equilibrates 
very fast compared to the Ising limit /3 = 0. 

The detailed dynamical studies in fact showed that the system belongs to the Model I 
dynamical class for any finite /? and a dynamic transition takes place exactly at /3 = 0. Certain 
subtleties in this model need to be mentioned: with respect to Model I, /3 = is the maximum 
noise and its inverse may be thought of an effective temperature. On the other hand, from 
the Ising model viewpoint, /? plays the role of noise. However it is not equivalent to thermal 
fluctuations which can affect the state of any spin. With /3, flipping of spins can still occur 
at the domain boundaries only. Hence, even with this noise, the equilibrium behaviour is not 
disturbed for any value of /? (even for /3 — )• oo which corresponds to Model I) while in contrast, 
any non-zero temperature can destroy the order of a one dimensional Ising model. 

2.2.3. Model with random cutoffs : Previously we have discussed the case when a uniform 
cutoff is introduced to the system. One can introduce randomness in the cutoff parameter p 
varying from to 1, (chosen randomly from a uniform distribution) and associated with each 
individual a different value of p |T7j. The randomness is quenched as the value of p assumed 
by any individual is fixed for all times. In this case, the variation of the magnetization, domain 
walls and persistence show power law scalings with exponents corresponding to Model I only 
for an initial range of time; however, calculation of z from the saturation times gives z = 1. 
The later time behaviour appears to deviate from the Model I behaviour and further studies are 
required to analyses the exact dynamical behaviour. 




Figure 7. Snapshots of the system in time for (a) random cutoff parameter p varying from 
to 1, and (b) uniform cutoff parameter p. These snapshots are for a L = 100 system. 

Even snapshots of the system (when compared to the uniform cutoff or infinite cutoff cases) 
do not help in understanding the phenomenon much (Fig [7|) . 

3. Mapping of the opinion dynamics model to reaction diffusion system 

The opinion dynamics model discussed so far are models where the dynamics is described in 
terms of the Ising spins that mimic the binary opinions an individual can have. In this model, 
a spin deep inside a domain does not flip. The dynamics is governed by the flipping of the spins 
only at the domain walls. The dynamics, in this respect, is reminiscent of the zero temperature 
Glauber dynamics of the kinetic Ising model. The motions of the domain walls can be viewed as 
the motions of the particles A with the reaction ^4 + ^ — t- 0. This means the particles are walkers 



and when two particles come on top of each other they are annihilated. The annihilation reaction 
ensures domain coalescence and coarsening. Unlike that in Glauber Ising model, the walkers A 
corresponding to Model I do not perform random walks. These walkers move ballistically towards 
their nearest neighbours. This bias, as we have seen before, gives rise to a new universality class 
than that of conventional reaction diffusion system [20] . 

We have also studied ^ + A — )• model with the particles A performing random walk with a 
bias e towards their nearest neighbors. We have taken e as the probability that a walker walks 
towards its nearest neighbour. Clearly, e = 0.5 corresponds to usual reaction diffusion system 
with the particles performing random walk. On the other hand, e = 1 is equivalent to our Model 
I as has been described above. We have studied the dynamics of the reaction diffusion system 
for different values of e in the range [0"^, 1.0]. 

In reaction diffusion systems, the growth of domains is given by the number of surviving 
walkers. Persistence P{t) in these systems is defined as the fraction of sites unvisited by any of 
the walkers A till time t. Figures [S] and O show the decay of persistence and fraction of walkers 
with time for different values of e > 0.5. 
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Figure 8. Decay of persistence with time for e = 0.7 and e = 0.55 

We find that for e > 0.5, the Model I behaviour is observed, namely: z ~ 1.0 and 9 ~ 0.235, 
with some possible correction to the scaling which becomes weaker as e is increased. For example, 
there is a logarithmic correction to scaling for the decay of the fraction of domain walls which 
takes the form t~^{l + a(e) log(t)) where a{e) — ?■ as e — >• 1. One can compare the above model 
with the cases discussed in section 2.2.3, where the introduction of stochastic dynamics also 
occurred with a bias towards the larger domain. In case of thermal disorder, the parameter 
comparable to e is /3. In the present case, the exact sizes of the domains do not matter (which is 
important for the case with f3) but the results are consistent in the sense that any bias towards 
the larger domain (or nearest walker) makes the system behave like Model I. 

In this model, we have also studied the e < 0.5 region where the opposite happens, the walker 
has a bias towards the further neighbour. Obviously domain annihilations take place very slowly 
now, even slower than l/log(i) and the dynamics continues for very long times. Consequently, 
the persistence probability no longer shows a power law variation now but falls exponentially to 
zero (Fig [10]). 
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Figure 9. Decay of number of walkers with time for e = 0.7 and e = 0.55. There is a logarithmic 
correction to scaling for the value of e = 0.55. The form of f(t) is t~^(l + Qlog(t)) with a = 3.92 
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Figure 10. Decay of persistence and number of walkers with time for e = 0.3 and e = 0.45. 
The form of g{t) is a(l/log(t)) where a is any constant. 



4. Summary and Discussions 

We have reviewed the dynamics of a recently proposed opinion dynamics model (Model I and its 
variations) [8]. Model I in one dimension belongs to a new dynamical universality class with novel 
dynamical features not encountered in any previous models of dynamic spin system or opinion 
dynamics. In the corresponding reaction diffusion system A + A — t- 0, we have introduced a 
probability e of random walkers A moving towards their nearest neighbors, e = 0.5 corresponds 
to the particles A performing unbiased random walks and the system belongs to the dynamical 
universality class of zero temperature Glauber Ising model. We find that for e > 0.5, the system 
still shows power law behavior of domain growth and persistence but with a universality class 
of that of Model I. For e < 0.5, the domain grows logarithmically and the persistence decays 
exponentially in time. 

We have discussed quite a few variations on Model I with quenched and annealed disorders 



and cutoffs in the interaction range. In presence of quenched disorders, hke the presence of 
few rigid spins which do not flip, the dynamics changes drastically. The dynamics ends up in 
heterogeneous equilibrium phases, which is fully disordered in the sense that no consensus can be 
reached as the order parameter goes to zero in the thermodynamic limit, with identical behaviour 
irrespective of the percentage of the rigid spins in the system. However, the dynamics shows 
extreme robustness against annealed disorder. We have also discussed the crossover phenomena 
in the dynamics in presence of cutoffs in the range over which a spin recognizes its neighboring 
domains. 

All the models discussed in the present article can easily be extended to higher dimensions 
and its universality class determined. Phase transitions occurring at non-extreme values of suit- 
ably defined parameters may also be expected in higher dimensions. 
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